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Abstract. We study the geodesies in a planar chessboard structure with two values 1 
and /3 > 1. The results for a fixed structure allow us to infer the properties of the Finsler 
metrics obtained, with an homogenization procedure, as limit of oscillating chessboard 
structures. 



1. Introduction 

In this paper we deal with optical paths in a dioptric material with parallel geometry 
and a chessboard structure on transversal planes. Further to a bidimensional reduction, 
we fix the optical features of the composite material in terms of its refractive index: given 
/3 > 1, let us define on [0, 2) x [0, 2) the function 



(1) apix,y) 



if ix,y) G [(0,1) X (1,2)]U[(1,2) x (0,1)], 
otherwise. 



and extend it by periodicity to a function defined on M? which we still denote by o/j. 

Hence, normalizing to 1 the speed of light in the vacuum, light travels in the dioptric 
material with a speed l/o/j. 

Since we are dealing with a system employing only refraction, Fermat's principle dictates 
that the optical paths between two points minimize the optical path length, which coincides 
with the time spent. Thus, in an homogeneous material, where the speed of light is 
constant, the optical path is a segment. Moreover, Fermat's principle leads to Snell's law 
of refraction, which completely describes the optical paths in layered materials (see [?] for 
a comprehensive introduction on the principles of optics). The explicit description of the 
optical paths in the chessboard structure becomes harder since, for example, no necessary 
condition prescribes their behaviour at corners. 

From a geometrical viewpoint, we are interested in the description of the geodesies in the 
Riemannian structure (M^,a^), that hereafter will be called standard chessboard structure. 
We shall refer to light squares or to dark squares, when, respectively, = 1 or = 

In the mathematical model, a virtual path emanated from the origin is a solution 
u: [0,T] ^ to the differential inclusion 



(2) 



iu'{t) G Gp{u{t)), 0<t<T, 

h(o) = 0, 

where the set-valued map (with nonconvex values) is defined by 

Gf3{x, y) = — ^ dBiiO), {x, y) G 
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(see e.g. [?] for an introduction to differential inclusions). 

Fermat's Principle states that a ray of light from O to ^ is a solution of the minimum 
time problem with target ^ 

(3) Tp{C) = inf{r > 0; 3 u{-) solution of (|), with u{T) = ^} . 

In [?] it is proved that the infimum in (^) is reached. Moreover, if n(-) is a solution of (P) 
satisfying u{T) = S^, then the optical length of the curve F = {u{t); < t < T} is 

Cp{T) = r ap{u{t))\u'{t)\dt = T. 
Jo 

Then optical paths are geodesic curves in the Riemannian structure (M^, 0/3). We underline 
again that the global minimization procedure considers only refracted rays, excluding all 
reflected rays, because a reflected ray is never a global geodesies. Hence our results have 
a reasonable physical meaning either if the number of interfaces traversed in the periodic 
media is not very high, or for /3 near to 1, since in both cases the reflected light can be 
neglected. Anyhow the results are intended as a depiction of curves of minimal length in 
the Riemannian structure (M^,a^). 

A starting point is the elementary observation that, for /? > 2, any geodesic joining two 
points in the light material (i.e. in the set {a^ = 1}) is a "light path", i.e. it never crosses 
the dark material. With little more work it is not difficult to prove the same conclusion if 
/3 > 1/2, provided that the endpoints of the geodesic have integer coordinates. Moreover, 
the value /3 = \/2 is an optimal threshold in this class of geodesies, since the diagonal of 
a dark square is a geodesic for every /3 < \f2. 

In this paper we obtain a perhaps surprising result, focusing our attention to geodesies 
joining two points with coordinates (2n + j, j), n, j E Z, that we call light vertices. For 
/3 > Y^3/2 we depict explicitly the geodesies joining the origin to a light vertex. As a 
consequence, we prove that the threshold value for the minimality of light paths is given 
by /3o G (^^3/2, \/2), which is exactly the value of /3 such that the optical path joining 
O with the light vertex (3, 1) has the same length (2 + of the optimal light paths. 
More precisely, we show that for /3 > /3g the geodesies are optimal light paths, whereas for 
a/3/2 < f3 < Pq the minimal curves are constructed concatenating the maximal number 
of translations of the optical path joining O with (3, 1), with segments either on the sides 
of the squares or on light diagonals. 

For 1 < /3 < a/3/2 the characterization of the geodesies seems to be a hard problem, 
for reasons that will be clarified in Section |5|. Anyhow, we are able to compute the optical 
paths joining the origin to a light vertex in a small cone {0 < K y < x}, where K = K{j3) 
is an odd positive integer, whose value diverges to +00 as (3 approaches 1. 

The knowledge of the minimal length of curves joining the origin to light vertices in the 
chessboard structure is enough to characterize the so-called homogenized model. Namely, 
the results described above give information about the optical paths in an inhomoge- 
neous dioptric material whose observed refractive index, at a mesoscopic level, is given 
by the chessboard structure, that is the observed refractive index at a given scale e > 
is ap{x,y) = ap{x/e.,y/e). We are interested in the behavior of the optical length of 
geodesies when e ^ 0. 
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We already know that, at the scale e, a virtual path emanated from the origin is a 
solution u: [0,T] ^ of the differential inclusion 

f^'(*) e dBi{0), 0<t<T, 

(4) <^ ' ' a'pHt)) 

[u{0) = 0, 

and Fermat's Principle states that a ray of light from O to ^ is a solution of the minimum 
time problem with target ^ 

(5) T^(0 = inf{T > 0; 3 n(-) solution of (|), with u{T) = ■ 

We are interested in the characterization of the limit, as e — > 0, of the minimum time 
problems (I)-®. 

The minimum time problems can be rephrased in terms of minimum problems of the 
Calculus of Variations. Let us denote by £^ the length functional in the chessboard 
structure corresponding to a|, that is 

£^(n) = j\^p{u{t))\u'{t)\dt, nG ^C7([0,1],m2)^ 

and by d^(0,^) the distance between O and ^ in such a structure, that is 

(6) d^(0, i) = inf : u G AC{[^, 1]; M^) s.t. u(0) = 0, = ^} . 

If u{-) is a solution of @ satisfying u{T) = ^, then T equals the optical length of the 
curve r = {u{t); < t < T} so that 

(7) T|(0 = 4(0,0, eeM^. 

The advantage of this formulation is that the asymptotic behavior of can be discussed 
in terms of F-convergence of the functionals (see e.g. [?] for an introduction to F- 
convergence) . In [?] it was shown that the sequence (£^) F-converges in 74C([0, 1], M^) 
(w.r.t. the topology) to the functional 

£^°'»(n)= ['<Pf,{u'{t))dt, uGACi[0,l],R''), 
Jo 

where — > [0, +cxd) is a convex, positively 1-homogeneous function, such that |^| < 
^/siO < for every ^ G M^. As a consequence, <Pp turns out to be a homogeneous 
Finsler metric in (see e.g. [?] for an introduction to Finsler geometry). 

In [?] it is shown that (Pp is not a Riemannian metric in for every /? > 1. In this 
paper we shall refine this result proving that the optical unit ball {'Pp < 1} is neither 
strictly convex nor differentiable (see Theorem 6.4 and Corollary |6.5| below). 

Since 'Pp is characterized by 

^^(e) = ^hm 4(o,e), 

then the limit of the minimum time problems (^-(^) is given by 

r^(0 = mf{r > 0; 3 u{-) such that <?/3(u'(t)) = 1, u{0) = 0, u{T) = ^} = 'Pp{C) . 

This is what we have called homogenized model related to the chessboard structure. 

Since is positively 1-homogeneous, it is completely determined by the geometry of the 
optical unit sphere {'Pjs = 1}, which is, in some sense, a generalized geometric wavefront 
with source point located at O. Moreover, due to elementary symmetry properties, it is 
enough to describe the set {^/j = 1} H {0 < y < x}. 
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Starting from our results on chessboard structures, we obtain that, if /3 > /3q, then the 
homogenized metric is 

(8) c|)^(x,y) = (^2- l)min{|x|,|y|} + max{|x|,|y|}, V(x,y)eM^ 

and the geometric wavefront {(^^ = 1} is the regular octagon inscribed in the unit circle, 
whose vertices lie on the coordinate axes and on the diagonals. This result generalizes 
the trivial remark concerning the case (5 > 2 (see [?]). On the other, we obtain that 
the octagonal geometry of the wavefront breaks for P = Pq, and the optical unit sphere 
becomes an irregular polygon with sixteen sides for a/3/2 < /3 < Pq (see Figure |T|). These 
are two of the main results of our paper; we refer to Theorem |6.3| below for their precise 
statement. 

For 1 < P < \/3/2 we will be able to compute the optical unit ball in the two small cones 
{-f^ |y| ^ l^;]} and {i^ |a;| < |y|}, where K = K{P) is the odd positive integer introduced 
above for the chessboard structure. More precisely, we shall show that in these regions 
the boundary {^p = 1} is piecewise flat with corners at the points (1,0), (0, 1), (—1,0), 
and (0, —1) (see Theorem |63^ ). As a consequence, the optical unit ball < 1} is neither 
strictly convex nor differentiable. 

We conclude this introduction with a warning on the physical interpretation of our 
results concerning the homogenized model. First of all, the geometrical optics approxi- 
mation is valid if the lengthscale e is much greater than the wavelength of light. If e is 
of the same order of magnitude of the wavelength, then we fall in the domain of photonic 
crystals optics, for which the full system of Maxwell equations must be considered. On 
the other hand, even in the range of geometric optics, our global minimization procedure 
considers only refracted rays, excluding all reflected rays. Since, at a macroscopic level, 
the number of interfaces traversed in the periodic media can be very high, the reflected 
light cannot in general be neglected. 



In the paper the following notation will be used. 

• |P, Q\. closed segment joining P, Q G 

. IP, Q] = [P, Ql \ {P}, [P, Q[ = IP, gi \ {Q] 

• |Pi, P2, . . . , Pn]: polygonal line joining the ordered set of points Pi, P2, . . . , Pn- 

• [t\ = max{A; G Z : k <t]. 

df 

• fs = -7^'- partial derivative of a function / with respect to s 

OS 

• S{A,B): Snell path joining A to B (defined in Section |2|) 

• >C/3(r) = Cp{T): length of F in the standard chessboard structure (optical length). 



2. Snell paths 
Let us consider the flat Riemannian structure (M^,a^), where 

- / \ f 1 if I xj is even, n n ^ 

ap{x.y) = i . . pen, p>i. 

yp if \x\ IS odd , 

Note that this structure corresponds to a composite medium whose structure is made 
by alternate vertical strips of light and dark material. 

It is well known that for every pair A = {xa^Va)-, B = {xbiVb) there exists a unique 
curve of minimal length (the geodesic curve) joining A to P, which is an affine path in 
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every vertical strip {[xj = fe}, G Z. Moreover, at every interface between two strips, the 
change of slope is governed by the Sneh's Law of refraction 

(9) sin 6*1 = /3 sin 6*2, 

where 6i and 02 are the angles of incidence with the interface from the light strip and the 
dark strip, respectively (see, e.g., [?, §3.2.2] or [?, §3.4]). 

In the sequel, this geodesic curve will be called the Snell path joining A and B and will 
be denoted by S{A,B). In order to fix the ideas, we shall always assume that xa < xb, 
and UA < Ub- We shall refer to the positive quantity xb — xa as the thickness of the Snell 
path. 

Let p, (7 > be, respectively, the thickness of the light and of the dark zone crossed by 
the path ^(A, i?), so that xb — xa = P + q, and let h be the vertical height h = i/b — Va- 
Since h = pianOi + gtan02 (see Figure ||), and (^ holds true, then a{p,q,h) = sin^i is 
implicitly determined in term of p, q and h by the constraint 



Clearly, a{p,q,h) is a strictly decreasing function w.r.t. p and q, while it is a strictly 
increasing function w.r.t. h. 

Finally, the optical length of the Snell path S{A, B), given by p/ cos 9i + I3q/ cos 02, can 
be expressed in terms of p, q and h taking again into account (^): 

(11) Lip,q,h):=^ P + ^''^ 



^2 



\/l - cr{p,q,hf - a{p,q,hf 

Lemma 2.1. Let a and L be defined by ( |7^ and (11) respectively. Then 

Lp{p, q, h) = Vl - 0-2 , Lq{p, q, h) = 
for every p, q > 0, and for every /i G M. 
Proof. Differentiating (|^) w.r.t. p, we get 

a pcTp P'^q^p _ 



+ 



(1-<t2)3/2 (^2 _ ^2)3/2 
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A. 


/ 













p(t) 

q(t) 



Hence 



Figure 2. The solid line corresponds to the Snell path S{A,B), while the 
dashed line is a minimal path joining A and B without crossing the dark 
squares. 



(1-^2)3/2 (/32_ ^2)3/2 

1 



a 



2 



2 



By an analogous computation one obtains the expression for Lq. □ 



Remark 2.2. By Lemma 2.1, it follows that, given the thickness t = p + q and the height 



G R of a Snell path we have that 
d 



—L{t - q,q,h) = -Vl - + J (3^ - > . 
dq ^ 

The geometrical meaning of this formula is clear: for Snell paths with fixed thickness, the 
more is the thickness of the dark material crossed, the more is the optical length of the 
path. 

3. The normalized length 

Let us consider now endowed with the standard chessboard structure. 

Definition 3.1. The n-th light diagonal, n G Z, is the straight line of equation 
y = X — 2n. A light vertex is a point having integer coordinates and belonging to a light 
diagonal. 

Definition 3.2. Given two points A and B in the same horizontal strip {y £ [r, r + 1]}, 
r G 7j, the Snell path joining ^ to i? (in the chessboard structure) is the geodesic S{A, B) 
in the corresponding parallel layer structure a{x + r,y). 

We are interested in the properties of the Snell paths starting from a light vertex A 
(say A = O, without loss of generality) and ending in a point B on the other side of the 
horizontal strip containing A (say i? = (x^, 1), > 0) (see Figure §). 

If < < 1, clearly S{0,B) = [0,i?] is the unique geodesic joining O to B. On 
the other hand, S{0,B) need not to be a geodesic when x^ > 1. Namely, we already 
know that for P large enough the optical length of S{0,B) is strictly greater than the 
optical length of the path obtained by a concatenation of horizontal segments on the lines 
X = or X = 1, with total length xs — 1, and a segment on a light diagonal, with length 
In this section we discuss the behavior of the difference Ci3{S{0, B)) — xb + 1 — \/2 
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Figure 3. Plot of f3 = 1.15. 



for > 1. To this aim, for t > and /? > 1, with some abuse of notation we define 
(T(t,/5) = a{p{t), q{t),l), where 



(12) 



lit) 



2 



if [t\ is odd, 
t — ^ if [t\ is even. 



and p{t) = t + 1 — q{t). Recah that a{t,fi) is determined by the constraint 



(13) 

As a consequence we have 



+ 



1 = 0. 



(14) 



V 



+ 



(5' 



< 0. 



t > 0, t 0N, 



(1-^2)3/2 (^2 _ ^2)3/2 

since pt(t), qt{t) > 0, and Pt{t)+qt{t) = 1 for every t > 0, t N . Thus the map t ^ a{t, fi) 
is strictly decreasing in [0,+oo), and satisfies (5"(0, /3) = l/\/2, linit_»+oo '5'(i, /3) = (see 
Figure ^). 

Moreover, as a straightforward consequence of the Implicit Function Theorem, we have 
that a^{t,(3) > for every t > 0. In particular we have 

(15) a{t,P)>a{t,l)= \ , Vt>0, V/3>1. 

V(l + t)"^ + 1 

Definition 3.3. We shall call normalized length the function 

Pit) , /32g(t) 



(16) 



l{t,P) :-- 



l-a{t,(3f 



+ 



/32-<T(t,/3)^ 



Notice that Z(t, f3) is the length of the Snell path joining the origin (0, 0) with the point 
(t + 1, 1) normalized by subtracting the minimal length of the paths joining the same two 
points without crossing the dark squares (see Figure ^). 

In order to simplify the notation we introduce the sets 

Il= [j {2k + l,2k + 2), lD=[j {2k, 2k + 1) . 

fceN fceN 

If t £ II then the last segment of the Snell path is in the interior of a light square, while, 
if t G Id, it is in the interior of a dark square. 
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The basic properties of the normaUzed length are collected in the following lemma. 

Lemma 3.4. The following properties hold. 

(i) kit,P) = ^l-a^{t,P)- 1 for every t £ h; 

(ii) kit, (3) = V/3^-<T^(t,/3)-l for every t G Id; 

(iii) l{-,l3) is strictly convex in any interval of II U Id; 

(iv) l{-,(3) is strictly monotone decreasing in any interval of II; 

(v) if (3 > \/3/2, then l{-,f3) is strictly monotone increasing in any interval of Id; 

(vi) if I < (3 < \/3/2, then there exists a unique to > 0, characterized by a{tQ,P) 
VP^ — 1, and such that 

it{t,p)<o, y te[o,to)niD, 

lt{t,(3)>0, y t£ito,+oo)nID■ 



Proof. The derivatives in (i) and (ii) follow from Lemma 2A_, upon observing that pt = 1 
and qt = in II, while pt = and qt = 1 in Id- Clearly (i) implies (iv) , while (ii) and the 
fact that < a\t,f3) < 1/2 imply (v) and (vi). 

(iii) follows from (i), (ii), and the fact that a{-,l3) is a decreasing function. □ 



In conclusion, since 1{0,P) = for every /3 > 1, by Lemma |3.4| we have that, for 
/3 ^ v^3/2, l{t,(3) > for t £ (0, 1) and the local minima of /(-,/?) are attained at t = 2k, 
A: G N, corresponding to the Snell paths ending in the light vertices (see Figures || and 

On the contrary, if (3 < ^3/2, a new local minimum for l{-, (3) may appear (see Figures ^ 
and 0). One may wonder if l{to,l3) is an absolute minimum for some /3. The following 
result shows that this is never the case. (We warn the reader that the proof is rather long 
and technical, and can be skipped in a first reading.) 



Theorem 3.5. Given 1 < (3 < ^3/2, let tQ > be as in Lemma \3.4 fvi). Then 

li2ko + 2,(3) <lito,f3), 
where ko = min{A; G N: to ^ 2A; + 2}. Moreover, the strict inequality holds if to ^ 2kQ + 2. 



Proof. If to G [2ko + 1, 2A;o + 2], then by Lemma |j(iv),(vi), l{2kQ + 2,(3) < l{t, (3) for every 
t G [0, 2/co + 2), and the result is straightforward. 

Let us now consider the case to G (0, 1), so that k^ = Q. 

Recalhng that a {to, (3) = x/p^ -1, we obtain that l{to, (3) = \/2 - + - 1 - ^2. 
On the other hand we have 

l{2,(3) = ^^ + ^^=-2-V2-as\^^^+ -1 




= 2^1 - al + ^(3^-al-2-V2-a^, 

where (T3 = (t{2,(3), and we have used the constraint (p^ satisfied by CJ3. Hence, denoting 
by 

(17) ^{a, /?) = 2^1^ +^fi2_^2_2-a-^2-(3^-^(3^-l, 

we have to show that l{2,f3) — l{to,(3) = 99(0-3, /?) < 0. 

One can easily check that (p{a,f3) is strictly monotone decreasing w.r.t. a in [0,1], so 
that, by we obtain (/j(cr3,/3) < (p{l/VlO,P) for every /3 > 1. In addition, ip{l/^/lO, f3) 
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is a strictly monotone increasing function w.r.t. /3, so that we get 

</.(cT3, /3) < ip{l/^, (3) < ip{l/Vw, ^2) = -A=-2 + yi?-^/2<0, 

for every /? G (1, ^3/2), which concludes the proof for to G (0, 1). 
Assume now that to £ (2^0) + 1) with > 1. 
Since we have 

/(2A:o + 2, /3) = /(to, (3) + / /t(t, /?) dt + / lt{t, p) dt 

J to J2kQ+l 
f-2fco+l / I \ /•2fco+2 



dt + jj^^ {Jl-a^(t,l3)-l^ dt, 



our aim is to prove that 

/•2feo+l , /•2fco+2 . 

(18) / ^(3^ -a^{t,f3)dt+ ^l-a^{t,f3)dt <2ko + 2-to. 



to J2fco+l 



We split the proof of (18) into three steps. 
Step 1. Setting 

/(A;,t,7) = 2A; + l + (/32-l)t + ^2 

(19) , / 1 + /l + _ i)e-2/{fc+i) 

/3v/(2A: + 2)2 + (t + l)2(/?2_i) - (k + I) log ( ^ 



and c = 1/(1 - a{2ko + l,/3)^), we show that if f{ko,to,c) > then (|T|) holds. 
Step 2. Setting 

(20) g{b) = l + 3b'- 46^ + (1 - e'^) ^^^^ , co = co(6) = 1 + ^{b' - 1), 

V^o 16 

then f{ko,to,c) > g{(3). 

Step 3. g{b) > for every b £ (1, y/3/2). 

Proof of Step 1 . Let us consider the functions 



Recalling (|l^, and taking into account that pt{t) = and qt{t) = 1 for t G [to,2A;o + 1) 
we obtain 



atCT ^ &^ 1 ^^y/F^ 

P^-^^' P f3h ^ /32(t + l) ' 

(1 _ <72)3/3 + (^2 _ ^2)3/3 

where in the last inequality we have used the fact that the function b ^ }? ji}? — o"^)^/^ is 
strictly monotone decreasing, and p + g = t + l. In conclusion we obtain that ^ satisfies 
the differential inequality 

(21) [^0,2^0 + 1), 

U(to) = l. 
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Similarly, recalling that ptit) = 1 and qt{t) = for t G (2/co + l, 2A;o + 2), we obtain 



ata 1 



(1 _ ^2)3/3 + (^2 _ ^2)3/3 

and, since p > ko + 1, we conclude that x satisfies the differential inequality 

(22) tG (2^0 + 1,2^0 + 2), 

[Xi2k0 + 1) = ^l-/?2+^2(2A;o + l). 

Solving the Cauchy problems associated to the differential inequalities (|2l|), (|2^), we get 

(23) <>(*)< , ' , ■ «€ [(0,2*0 + 11, 



/32 V /3V (i + l)2 
(24) x{t) < I ^ t€[2ko + l,2ko + 2]. 



,l-/32 + V? 
As a consequence of these estimates we obtain 

- a^it, I3)dt<(3 ^^{2k^ + 2)2 + (to + l^ifi^ - 1) - /3(to + 1) 

/.2fco+2 / / 1 + Jl + (c - l)e-2/{'=o+i) 

/ Jl-a2(t,/5)dt< l + (A:o + l)log ^ ^— 



where c = — 0^ + ip{2kQ + 1)2), which concludes the proof of the Step 1. 
Proof of Step 2. From ([23|) and the fact that to > 2ko, we have that 



+ 1) < .,„, jr-"^:;-^ < + ^) 



v/(2fc„ + 2)» + (/3^-l)((o + l)'' v/(2fc„ + 2)» + (/3^-l)(2k„ + iy'' 

so that 



1 + (/?' - 1) 

> 



2fco + 2 



1 _ /32 + ^(2A;o + 1)2 - /2A:o + l 



2 



Moreover, it can be easily checked that the function 



1 + Jl + (7 - l)e-2/('=o+i) 
7 ^ log I -^=^^ I , 7 > 1 



is strictly monotone decreasing, while the function 



t^ (/?'-l)t-/?V(2fco + 2)2 + (t + l)2(/?2_i), t e (2A:o,2A:o + l) 
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is strictly monotone increasing. Hence we have that 

f{ko, to, c) > f{ko, 2ko, com = 2A:o/3' - P ,J {2ko + 2)^ + {2ko + 1^ {P^ - I) 

, 1 + 7l + (co-l)e-2/(fco+i)^ 

+ l + /32-(A;o + l)log ^ — — 

Vco + 1 



where cq = co(/3) is defined as in 
In addition, the functions 



k ^ 2k0^ - (5^ {2k + 2)2 + {2k + l)2(/32 _ 1) 
k ^ {k + 1) log 



1 + + (co - l)e-2/(fe+i) 



are strictly monotone increasing for k>l. Hence we get 



f{koM, c) > 1 + 3/?2 - f3JlQ + 9(/32 - 1) - 2 log 



'1 + Vl + (co-l)e-^' 



^ + 1 
Finally 



^ / 1 + y/l + (cQ - l)e-i \ + (cQ - l)e-i - ^ 



1 + (1 - e-i)i^ - 1 1 1 - /?^ 

CO <7T(l-e ) 



so that the proof of Step 2 is complete. 
Proof of Step 3. We have that 



g'{b) = 6b 



7 + 18^2 36(1 - e~^)6 
V7 + 9F^ (7 + 962)3/2 ' 



and, for l<b< 



A l^fi/, _ 7 + 18fc^ \ _ ^ 7 + 662 _ 36 7 + 662 

db[^^~V7 + W) ~ (7 + 962)3/2 - ^ - ^7 + 962 7T962 

>6-9^1^>0, 
a/41 16 ' 



and 

d r 6 \ _ 1862 -7 



16- r7T^<0, 6>1. 



d6 V (7 + 962)3/2; (7 + 952)3/2 

Then we get 

5'(6)>6-| + ^(l-e-i)>0, Vl<6<^. 

Hence ^'(6) > g{l) = for all 1 < 6 < ^3/2, and Step 3 is proved. □ 

Now we focus our attention to the study of the sequence l{2k,l3), € N. Given k £N, 
we set S{k, f3) = l{2k + 2,13)- l{2k, (5), that is 

(25) i(,,;3) = 4±L + '*' " 



■/r^ v^?^^ v/r^^ v^p^^ 
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where r = r(/c, /3) := a{2k + 2, (3) and fi = a{k, (3) := a{2k, (3) are imphcitly defined by 

(k + 2)T (k + l)T 

(27) (^ + ^)^ + - 1 

By the monotonicity of the function o" (•,/?) (see inequahty ([l^ ) it fohows that r < ex. 

Remark 3.6. While the sign of l{t,l3) gives a comparison between the optical lengths of 
the SneU path S{0, {t + 1, 1)) and the "light path" [O, (1, 1), (1 + t, 1)], the sign of S{k, (3) 
gives a comparison between the optical lengths of 5'(0, (2A;+3, 1)) and of 5'(0, (2/c+l, 1))U 
|(2/c + l,l),(2A: + 3,l)l. 

Remark 3.7. Given /3 > 1, consider the function /: [0,+oo) — > M, affine on each interval 
[2k,2k + 2] and such that l{2k) = l{2k,P), ken. Then the derivative of l{t), for t G 
{2k,2k + 2), is given by 6{k,(3)/2. 

Since 6{k, (3) = Cp{S{0, {2k + 3, 1))) - Cp{S{0, {2k + 1, 1)) - 2, and it is clear that 
C0{S{O,{2k + 3,1))) - C(3{S{0,{2k + 1, 1)) ~ /3 + 1 for A; ^ +oo, one expects that 
5{k,P) ~ /3 — 1. A more precise result is the following. 

Theorem 3.8. Let (3 > 1 he fixed. Then {5{k,f3))i. is a strictly monotone increasing 
sequence and 

(28) ,(,,^, = (^-l,-_^L_^J,+0(J3), .^+00. 

Proof. We can define r, o", and 6 respectively through (^), (^) and ( p5| ) as smooth 
functions of E M, A: > 0. Differentiating 5 w.r.t. we get 

A ri- m = 1 , (fe + 2)TfcT , (fc + l)/3Vfcr 

0UK,p; + _ ^2)3/2 + y;52—^ + (/52-r2)3/2 

1 (A; + l)(Tfc(T /32 A:/3VfcCT 



(l-a2)3/2 ^/32 _ ^2 (^2 _ ^2)3/2- 

On the other hand, differentiating the constraints ( p6|) and (27) we obtain the identities 
(29) (fc + 2)T, ^ (fc + l)/?Vfc _ r T 



(1-t2)3/2 (/32_ ^2)3/2 ' 

(fc + l)crfc _^ kp'^ak a a 



(1-^2)3/2 (;52_ ^2)3/2 



2-^2 



Hence, being r < cr, we get 

6k{k,(3) = Vl-T^-\/l-<T^ + - r2 - y^/32 - ^2 > 0. 

In order to determine the behavior of 6{k,l3) for A; large, notice that, setting e = 1/A; 
and (7{e) = cr(l/e, /3), (p7|) becomes 

{l + e)a a _ 

^ V/32 - cF2 ~ ^ ' 
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that is a is implicitly defined by 



+ 



1 



One has /(O) = 0, /'(O) = /"(O) = 2^, so that 

SF(0) = 0, 5'(0) = -^, 5"(0) = -2- 



and hence 



(31) 5(.) = ^._^.^ + 0(a 

that is 

and 

.(.,«=a(. + l,/3)=;f(^)=.(i-J, + 0(^ 



/? + 1 A; 1 (/3 + l)2y fc2 ' y^s 

Finally we have 



/3) = (A; + 2) (^1 + -r^ j + + 1) (^1 + _ j - (fc + 1) (^1 + -a' 

completing the proof. □ 
Definition 3.9. Given /3 > 1, we shall denote by kc{P) the integer number defined by 
(34) kciP)=mm{keN: 6{k,P)>0}. 

By the very definition, we have that l{2kc{P), (3) < l{2k,(3) for every /c G N (see also 



Remark 3^). Moreover, by Lemma |3^ and Theorem |3.5| , the absolute minimum of 
is attained at a point t = 2k, k £N. Hence l{2kc{P), (i.e., the normalized length of the 
Snell path joining the origin with the right-top vertex of the {2kc + l)-th square) minimizes 
the normalized length among all the paths remaining in a single horizontal strip. 

Now we want to study kc{P), /? > 1. As a preliminary step we investigate the behavior 
of 5{k, •) for a given k. 

Lemma 3.10. Let k €zN be fixed. Then the function jS ^ 5{k,P) is strictly increasing in 
[1, +oo) . Moreover 6{k, 1) < and 6{k, \/2) > 0. 
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Proof. Differentiating 5 w.r.t. (3, we get 

6f3{k, 13) = 1^ _ J^|3/2 + _ ^2)3/2 + + 1) (-^2^_ ^2^)3/2 

^^^^ _ {k + _ /?^-2a2 _ 

(1-(t2)3/2 ^ (^2 _ ^2)3/2 (^2 _ ^2)3/2 " 

Differentiating (p6|) and (|27|) w.r.t. /?, we obtain 
(36) £±^ + ,t + l, '^I!!^ = (, + 1). 



(1 _ ^2)3/2 ' ^ ^ (^2 _ ^2)3/2 ^ ^ (^2 _ ^2)3/2 

(37) P^ + k-f^ = k- 



(1 - (72)3/2 (/32 _ ^2)3/2 (^2 _ ^2)3/2 " 

Substituting (^) and ( |37|) into (|35|), we conclude that 

We have to show that 
(38) iiiiE - > 0, 

For every k G M, let us denote by s{k) the unique function implicitly defined by 
(391 (^ + ^)" + - 1 



so that s{k) = a, s{k + 1) = r. Since inequality (pq) clearly holds true for /c = 0, it is 
enough to show that 

d I K \ /32 — §2 _|_ f^gg 



Differentiating (|39|) w.r.t. «; (see also (^9|) ), we get 



1 1 

+ 



Sk(«^) = -" 



K + 1 



(1-s2)3/2 (^2_g2)3/2 

Moreover, using again (|39[), we have 



so that 



1 1 _ 1 1 1 

Vl - s2 \//32 - s2 kVI - s2 

1 



+ 



- + KSSk = /32 _ ^2 _ 



2_g2 



K+1 _^ 



(l-s2)3/2 (^2_52)3/2 

K/32 



> ( ^/?' - - ) >o, 
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where the last inequahty can be easily checked recalling that < s < 1, while /?, k > 1. 
Hence we conclude that inequality ( ^8|) holds true, which implies that the function 5{k,(3) 
is strictly increasing w.r.t. (3. 

Taking into account that, by (p7|), (j(0, \/2) = \/2/2, we easily get that 5(0, a/2) > 0, so 
that, by Theorem U, 6{k, V2) > 0. 

In order to prove that 5{k, 1) < 0, we note that, from (|26|) and (^), we get 



t(/c,1) = -==i==, a{k,l) ^ 



so that 



5{k, 1) = V 1 + + 3)2 - ^1 + (2A; + 1)2 - 2 < , 
concluding the proof. □ 



Remark 3.11. As an easy consequence of Lemma 3.10 we obtain that kc{(3) is a nonin- 
creasing function of /?. 



Thanks to Lemma 3.10, the following definition makes sense. 

Definition 3.12. For every A; G N we shall denote by (3^ the unique number in {l,\/2) 
such that 6{k,(3^) = 0. 

By Theorem |3.8| we have that 

5(j,/3^)<0, \fj<k, 5{j,Pt)>0, yj>k, 

hence 

(40) l{2k, (31) = l{2k + 2, (31) < l{2j, PI) V j E N \ {k, k + l]. 
In particular we have that 

(41) l{2k + 2, (31) = l{2k, (31) < l{0, (31) = 0, 
where we have taken into account that l{0,(3) = 0, for every /? > 1. 

Lemma 3.13. The sequence (/?^)^gi^ is strictly decreasing and limfc^+oo/?fc = 1- 



Proof. Given /c G N, by the monotonicity of 5 w.r.t. k stated in Theorem 3.8 we have 

6{k,p^,_,,)<6{k + l,Pl^,) = 0, 

where the last equality follows from the very definition of Again, the definition of 

and the monotonicity of S w.r.t. (3 stated in Lemma 3.10 imply that P^^i < (3^- 
In order to prove the last part of the thesis, given A; G N, let us define the functions 



f'is, (3) ■■= ^Ji4 + -ji^ -2k-V2, 



^(3"^ - s2 

(42) {k + l)s ks _ 

h^{s,(3) := f^{s,(3)- sg^{s,(3) = {k + 1)^1 - s'^ + k^ (3"^ - - 2k - V2 + s , 

where s G [0,1) and /? > 1. Since a{2k,(3) is the unique solution of g^{s,(3) = 0, we 
have that {a {2k, (3), (3) = f''{a{2k, (3), (3) = l{2k,(3). Moreover /i^(s,/3) = -g^{s,0) and 
(7g (s, (3) > 0, so that h^{-,(3) is a strictly concave function in [0,1] which attains its absolute 
maximum at s = a{2k,(3). Hence 

l{2k,l3) = h''{a{2k,(3),p) = max h''{s,P) > /i^(0,/3) = {P - l)k + 1 - V2 , V/? > 1 . 

se[o,i) 
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Then, from (^), we have 

> l{2k + 2, p^) > {(31 - l){k + 1) + 1 - 

so that 

- 1 



(43) 

and the conclusion follows 



k + l 



□ 



The first values (up to the fifth digit) are listed in the following table. 



k 





1 


2 


3 


4 


5 


6 


7 




1.24084 


1.06413 


1.02820 


1.01577 


1.01006 


1.00698 


1.00512 


1.00392 



By Lemma |J(v), for (5 > y/3/2 the Sneh path S{0, {t, 1)), with t E (0, 1) is never a 
geodesic. This property will be crucial for the results in Section |5[ The following result 
gives the position of ^3/2 ~ 1.22474 w.r.t. the critical values /?^. 

Lemma 3.14. Pf < ^/3/2 < /3g. 

Proof. From (^) we have 



In order to prove the inequality \/<i/2 < (3q, let g^, h} be the functions defined in (42) 
for /c = 1, and let ai G (0, 1) be the unique solution of 5^(s, /Jq) = 0. Since 

= 5(0,/3§) = /(2,/3S) = /n(7i,/3S), 
we have that /i^((Ti,/3q) = 0. Moreover 



> 0. 



Hence the inequality \/3/2 < Pq can be obtained showing that 

-s^-2-V2 + s<0, 



(44) 



s2 + 



The inequality (44) easily follows observing that 



2\/l - s2 + -s < -VW, 
3-3 



3 „ 1 /5 

+ -s < \ -, 

2 3 - V 3 



Vs G (0,1). 



Vs G (0, 1) , 



so that 



<^Vw + ^-2~V2<0, VsG(0,l) 



which completes the proof. □ 

We summarize the previous analysis in the following result, which is depicted in Fig- 
ures 

Corollary 3.15. For every /? > 1 

(45) mml{s,p) =mml{2k,P) <l{t,P) ^ t £ [j {2k,2k + 2). 

Moreover the following hold. 
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Figure 4. Plot of (3 = 1.26 {(3 > 





Figure 5. Plot of lit, (3), (3 = 1.23 (73/2 < (3 < (3^) 





Figure 6. Plot of (3 = 1.2 < /3 < ^/3/2) 



(i) A;c(/3) = for every (3 > I3q, and kc{f3) = k + 1 for every /? G {(3^+1^ ^k\' ^ ^ 

(ii) ifP> (3^, then = 1(0,(3) < l{t,(3) for every t > 0; 

(iii) ifPG ken, then l{2k + 2,P) < l{t,(3) for all t £ [0, +oo) \ {2A: + 2} ; 

(iv) l{2k, (31) = l{2k + 2, (31) < l{t, (31) for every ken, and t e [0, +oo) \ {2k, 2k + 2}; 

(v) if (3> v/372, then l{t, p) > for t e (0, 1] . 

(vi) for every (5 > I, l{2kc{(3) , (3) < l{t,(3) for every t > 2kc{(3). 



Proof. Formula (^5|) summarizes the results in Lemma |3.4| and Theorem 
Let us now prove (i) . By Lemma 3.10| and Theorem |3.§| we have 



(46) 5{j,(3)<5{k,P)<5{k,Pl] 

(47) 5{3, (3) > 5{k + l,(3)> 6{k + 1, (3',+,) = , 



0, yj en, j <k (3<(3'k, 
VjGN, j>k + l,P>(3',^,, 
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Figure 7. Plot of l{t,(3), (5 = 1.009 < /? < /?|) 

and we have the strict inequahty in (^) if /3 7^ /?^. Hence, recalhng the definition of kc{P') 
given in (p^), we conclude that (i) holds. 

Moreover, since for every n,mGN,n<m, we have 

m—l 

l{2m,(3)=l{2n,(3)+Y. 5{j,(5), 

j=n 

as a consequence of (|^) and (p7|), for every /c G N we get 

(48) = l{0, 13) < l{2j, 13), Vi G N, j / 0, /5 > /3§ 

(49) /(2fe + 2,/3) < K2i,/3), Vj G N, j / A: + 1, /3 G (/J^+i,/?^ 

(50) l{2k,(3l) = l{2k + 2,(3l)<l{2j,(3l), Vj G N, j / A;, + 1, . 

Hence (ii) , (iii) , and (iv) follow from (^) . Finally (v) follows from Lemma |3.4| (v) , whereas 
(vi) is a direct consequence of (i)-(iv). □ 

Up to now we have described the behavior of the normalized length of Snell paths 
starting from a light vertex and remaining in a horizontal strip. The following result deals 
with the normalized length of any Snell path remaining in a horizontal strip. 

Lemma 3.16. Given x G M, r > 0, r G let S{A, B) he the Snell path joining A = (x, r) 
to B = (x + T, r + 1) . Then 

(51) C^iSiA, B))-T + l-V2> li2k,iP),(3) . 
Moreover 



i) if P ^ Pf. for every /c G N, then the equality in (51) holds if and only if r 
2kc{(3) + 1, and x = 2n, n G Z; 



ii) if [3 = for some /c G N, then the equality in ( \5]\ ) holds if and only if t E 
{2kc{P) - 1, 2kc{P) + 1}, and x = 2n, n£Z. 



Proof. By Remark |2.2| we have that 

(52) Cf,{S{A, B)) > Cp{S{0, C)), C = (r, 1) , 

since ^(O, C) crosses a quantity of dark material not greater then the one crossed by any 
other Snell path with thickness r . On the other hand, if r G [0, 1), then 



Cfs{S{0, C)) = Vl + r2 >r-l + V2, 

so that, by (52), 

C/siSiA, B))-T + l-V2>0> l{2kc{l3) + 2, /3) . 
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Moreover, we stress that the equahty in ( |5l|) never occurs when r G [0, 1). If r > 1, then 

Cp{S{d,C))=l{T,P)+T + V2, 

so that 

(53) Cp{S{A, B))-T + l-^/2> 1{t, f3) > l{2k,{P),(3) . 

It remains to discuss the occurrence of the equahty in ( |53| ) for r > 1. 

If /3 7^ f3^, G N, then has its strict absolute minimum point at t = 2kc{P), so 

that 1{t -!,/?) = I{2kc{[3),l3) if and only if r = 2kc{l3) + 1. 

Moreover, for r = 2kc{(3) + 1, the equality C(i{S{A,B)) = £/3(S'(0,C)) holds if and 
only if X = 2n, n G Z. Namely, if denotes the thickness of light material crossed by 
S{A, B), then px < kdP) + 1 = p{2kc{P) + 1) (since the total thickness is 2kc{(3) + 1), and 
Px = kciP) + 1 if and only if A is a light vertex (i.e. x = 2n, n £ Z). 

li P = for some A: G N, then the conclusion in (ii) follows from (|5^ ) and the fact that 
the absolute minimum of /(•, /3^) is attained both for t = 2kc{P) — 1, and t = 2kc{(3) + l. □ 

We conclude this section by stating some properties, which will be useful in the second 
part of Section ^, of the following generalization of the normalized length introduced in 
dH). Let q{t), t > 0, be the function defined in (|l|). For < /i < 1, let p{t, h) = t+h-q{t). 
Given /? > 1, let a{t,P, h) be the unique solution of the implicit equation 

(54) 4^ + ^^21^-^ = 0, 
and let us define the function 

(55) [(t, /3, h) = ^ ^(^'^^ + ^ ^''^(^^ -t-hV2. 

^l-a{t,P,hf ^P^-d{t,P,hf 

The function / is the normalized length of a Snell path starting from the point {—h, —h) 
and ending in (t,0). It is straightforward that a{t,P,l) and l{t,P,l) coincide with the 
functions a{t, /?) and l{t, /3) defined in ([l3|) and ([l6|) respectively. The function t l{t, P, h) 
has the same qualitative properties of 1{-,P, 1) studied at the beginning of this Section. 
More precisely the derivative 

'^l-a^{t,P,h)-l, [ftGiL, 



(56) ltit,P,h) 



^p'' -a^t,P,h) -1, ate Id, 



is a monotone increasing function both in the set II and in the set Id (see Lemma 2.1). 
Concerning the derivative Ih. w.r.t. h, we have the following result. 

Lemma 3.17. The function h l{t,P, h) is strictly convex and monotone decreasing in 
(0, 1] for every t > 0, P > 1, and 



(57) !h{t, P, h) = ^l- d^{t, p, h) + <T(t, /?, h) - V2. 

Proof. Since Phit, h) = 1 we have that 



Differentiating ( |54[) w.r.t. /i, we get 

^ ^ (1 - ^2)3/2 ^ (/J2 _ ^2)3/2 
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Figure 8. 



Substituting ( |59| ) in (^), we obtain (|^). It is straightforward to check that the function 
a 1-^ Vl — <5"^ + a is strictly increasing in (0,l/\/2). Therefore, since t > imphes 
a < 1/V2, it fohows that lhit,(3, h) < 0, for h G (0, 1). Moreover, by @) it follows that 
fj/j > 0, so that the function h i— > a{t, (3, h) is strictly increasing in (0, 1) for every t > 
and (3 > 1. Hence, is strictly increasing for h £ (0,1), which implies that l{t,(3,h) is 
strictly convex w.r.t. h. □ 



4. General properties of the geodesics in the chessboard structure 

Up to now we have investigated the properties of a geodesic joining two points on the 
sides of one horizontal strip in the chessboard structure. In this section we will study the 
properties of a geodesic starting from the origin O, crossing an arbitrary large number of 
horizontal strips, and ending in a light vertex (2n + j, j), n, j G N. If n = or j = , then 
the unique geodesic from the origin to the point (2n + j,j) is the segment joining the two 
points. Hence we shall further assume that j,n > 1. 

Throughout this section we shall assume that 

(60) r is a geodesic from the origin to the point (2n + n,j G N, n,j > 1. 
The basic properties of T are listed in the following two propositions. 



Proposition 4.1. Let F be as in (6C ). Then the following properties hold. 

(i) Let H be a closed half plane such that dH is either the line x = k, or y = k, or 
a light diagonal D^, for some k £ Z. Let F C T 6e a path, with endpoints A, B, 
such that f <Z H and A, B e dH. Then f = {A, Bj. 

(ii) Let A = (xatUa) G T and let , r+ C F 6e the two paths joining O to A and A 
to (2n + j, j) respectively. Then the following bounds hold: 

- if A £ Z X Z then F~ C [0,xa] x [0, y^] anrf F+ C [xA,2n + j] x [yA,j]; 

— if in addition A is a light vertex, then 

T- C {(x, y) G : < y < yA, y < x < y - yA + xa}, 
r+ C {(x, y) G : yA<y <j, y - VA + xa < x < y + 2n} 
(see Figure 

(iii) Let Q be the interior of a light or a dark square, and assume that FnQ 7^ 0. Then 
r n Q is a segment. As a consequence, F = U^]^|Pj_i, Pj], where Pq = (0,0), 
P/v = (2n + Pi 7^ Pk for i ^ k, and Pi belongs to the boundary of a square 
for every i = 0, . . . , N . 
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(iv) Let 9i denote the oriented angle between the horizontal axis and |Pj_i,i-*j], i = 
1, . . . ,N. Then 9i £ [0, 7r/2]. Moreover, if either Pi-i or Pi is a light vertex, then 
e^ G [0,7r/4]. 

Proof. Property (i) is a straightforward consequence of the local minimality of F and the 
fact that a segment \A, B\ contained in the lines x = k, or y = k, or D^, for some k £ Z 
is the unique geodesic from A to B. 

In order to prove (ii), we first observe that 

r C W = {{x, y) G : 0<y<j, y<x<2n + y}. 

Namely, if this is not the case, there exists an open half plane H such that if H W = 0, 
H CiT / 0, and dH is one of the lines y = 0, y = j, Dq, Dn, a contradiction with (i). 
Then, if A G T C W has the stated requirements, then the bounds in (ii) can be obtained 
reasoning as above with half planes with boundary given by a line of the type x = xa, or 
y = yA, or D^^^_y^y2- 

Property (iii) is a necessary condition for minimality, see, e.g., [?, Section IV]. 

In order to prove (iv), we notice that, by (ii), 9i G [0, 7r/2] whenever either Pi-i or 
Pi is in Z X Z, and 9i G [0,7r/4] if either Pi-i or Pi is a light vertex. In particular 6i, 
On G [0,7r/4]. Hence we have only to show that if G [0, 7r/2] and Pi-i ^'Lx'L, then 
Oi G [0, 7r/2]. This follows from the fact that Pi-i is in the interior of a side of a square, 
so that the Snell's law (P) holds. □ 

Remark 4.2. Since = 1 on the boundary of the squares, then Proposition |4.1| (iii) can 
be improved observing that the intersection of F with the closure of a light square is a 
segment. 

In what follows we will be interested in the intersections F n Dj., k = 1, . . . ,n. 



Proposition 4.3. Let F be as in (6C), and let Pi, i = 0,...,N be as in Proposition 
\^.\ (iii). Then the following properties hold. 

(i) For every k = 0, . . . , n there exist < r]k < Ck ^ j such that F n = {Ak, B^J, 
Ak = {2k + %,r/fc), Bk = {2k + CkXk)- Moreover Aq = (0,0) = Pq, Bn = 
{2n + j, j) = Pn, and Ck < Vk+i for every fc = 0, . . . , n - 1. 

(ii) If Ak 7^ Bk then rjk and Cfc are integers (that is Ak and Bk are light vertices). 

(iii) Given k = 1, . . . ,n, let i = 1, . . . N be such that Ak G J-Pj-i, PiJ. Then < Oi < 
7r/4. If in addition Ak is not a light vertex, then 6i ^ 0. The same properties hold 
ifBk G [Pi-i^Pil 



Proof. It is clear that F n Dk 7^ for every k = 1, . . . ,n. Moreover, by Proposition 4.1(i) 



and Remark |4.2| , F n Dk is either a single point, or a segment joining two light vertices. 
The inequality Qk ^ Vk+i, k = 0, . . . ,n — 1 and property (ii) then follow from Proposition 
Hiv). 

Let us now prove that (iii) holds. By Proposition ^j|(iv ), we know that Oi G [0,7r/2]. 
Moreover, if Ak is a light vertex, then by Proposition [4.l| (ii) with A = Ak we get Oi G 
[0, 7r/4],. Finally, it has to be Oi 7^ tt/4 otherwise ^ F n Dk, in contradiction 

with the definition of Ak- 

Assume now that Ak belongs to the interior of a light square. Then, by (ii), F n Dk = 
{Ak}, so that Oi < 7r/4. Finally Oi ^ 0, otherwise F has to be an horizontal segment, due 
to Snell's Law (|). □ 

Definition 4.4. Given A; = 0, . . . , n, we say that F cuts the light diagonal Dk if the points 



Ak and Bk, defined in Proposition 4.3, coincide and belong to the interior of a light square. 
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Figure 9. Construction of T, kc{P) = 1. 



For every r = 1, . . . , j, we shall denote by F^ the curve 
(61) Tr = Tn{r-l<y <r}. 

By Proposition |4.1| (i), the intersection A := H {y = r — 1} is a single point, as well as 
ioT B := Tr n {y = r} . The curve F^ is a Snell path joining the two points A and B, and 
lying in a single horizontal strip. 

Remark 4.5. In what follows we shall assume, without loss of generality, that Fi is a 
Snell path starting from the origin. Namely, if this is not the case, Fi = S{Co,Ci) 
where Cq = (xo,0) and Ci = < xq < xi. Let us denote by pi, qi respectively 

the thickness of the light zone and of the dark zone crossed by Fi , and by p2, q2 the 
analogous quantities for the Snell path S{0,C2), C2 = {xi — X2,l). We have pi + qi = 
P2 + Q2, and p2 < pi, so that, by Remark £/3(Fi) > Cp{S{0,C2)). Hence the curve 
5(0, C2) U [C2, Ci] U (F n {y > 1}) is a geodesic. 

The following result is another fairly general property of the geodesies based on the 
behavior of the function l{t,(3) studied in the previous section. 

Proposition 4.6. Let P > 1 be given, and let kdfi) be as in Definition (see also 



Corollary 3.15(i)). For every r = 1, . . . ,j the curve F^ defined in intersects at most 
kc{P) light diagonals. 

Proof. Set, as above, A = Tr Ci {y = r — 1} = (x^, t — 1) and B = Tr r\ {y = r} = {xB,r) 
respectively the starting and the ending point of the Snell path F^. 

Assume by contradiction that F^ intersects more than kc{P) light diagonals, so that 
to := XB - XA-l > 2kc{/3). 

Let A; E N be the smallest integer such that 2/c + r — 1 > xa, and let m G N be the 
largest integer such that 2m + r < xb- Set Aq := (2n + r — 1, r — 1), Bq := (2m + r, r). It 
is clear that Aq and Bq lie respectively on the first and the last light diagonal intersected 
by Tr (see Figure]^, so that, by assumption, m — n> kc{(3). 

Let us denote by = 2k + r — 1 — xa and Ab = xb — 2m — r. It is not restrictive to 
assume that < A^- The points Ai = (a,r — 1) and Bi = (a',r) defined by 

(Tl) Ai = Aq, Bi = {2m + r + AA + AB,r), if either < Aa < A^ < 1, or < Aa < 
1, 1< Ab < 2, Aa + Ab< 2; 

(T2) Ai = (2k + r - 1 - Aa - AB,r - 1), Bi = Bq, if either 1 < A^ < A^ < 2, or 
< Aa < 1, 1< As < 2, Aa + As > 2, 
are such that a' — a — 1 = to > 2m — 2k, and either Ai or Bi is a light vertex. Moreover, 
from a direct inspection we can check that the thickness of the dark zone from Ai to Bi 
is not greater that the one from A to 5, so that by Remark we conclude that 

(62) Cf3iSiA,,Bi))<CpiSiA,B)). 

Let us consider the case (Tl), so that Ai = Aq, and let B' = {2k + 2kc{(3) + r,r). 
The assumption m — k > kc{P) implies that B' lies on the left of Bq (possibly the two 
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points coincide). Let us consider the following new paths: T' = S{Ai,B') U |i?',-Bi], and 
Fi = S{Ai,Bi). Since 2kc{f3) < to, from Corollary ^.15| (vi) we deduce that 

Cp{T') = Ii2kc{(3),f3) +to + V2< /(to, P)+to + V2 = C^iTi) . 

Finally, setting f = U S{Ai,B') U lB',Bj and noticing that the segments lA,Aij 

and l-B, BiJ have the same length, we have that 

£;3(f) = £/3(F') </:/3(Fi) <£;3(F,), 

where the last inequality follows from (|6^). 

The analysis of the case (T2) can be carried out in a similar way, with obvious modifi- 
cations. □ 



5. Geodesics of the chessboard structure (/? > y/3/2) 

In this section we shall restrict our analysis to the case P > y^3/2, and we shall provide a 
complete description of the geodesics joining two light vertices in the chessboard structure. 
The case /3 < ^3/2 seems to be harder to characterize, as we shall show by an example 
(see Example |5.6| below). 

The next theorem states that, for /3 > (3^, any geodesic F joining the origin with the 
point (2n + is a finite union of segments, connecting light vertices, and lying on light 
diagonals or on horizontal lines. 



Theorem 5.1. Let T be a geodesic as in ^duj). If (3 > /Sg, then the points A^, B^ are 
light vertices for every k = 0, . . . ,n, and {Bk-i, Aj.} is an horizontal segment for every 
k = 1, . . . ,n. 

Proof. By Corollary p.l5| (i) we have that kc{(3) = 0, so that, by Proposition f4.6| , F never 
cuts a light diagonal. Hence the points A^, B^ defined in Proposition |4.3| (i) are vertices 
of light squares, for every k = 0, . . . ,n. Given r = 1, . . . , n, let Bj—i be the exit point 
from Dr-i and Ar be the access point to Dr, and let i, i' = 1, . . . N he such that Br-i G 
lPi^i,Pil, and Ar £ We have that 9i = 6*^/ = 0, that is both [Pi_i,Pi] 

and |Pj'_i,Pj'] lie on the horizontal sides of the squares. Indeed, by Proposition ^l^(iii), 
6i',9i £ [0, 7r/4), and, if either 9i or 9ii belong to (0,7r/4) , then F should contain a Snell 
path starting from a vertex of a light square and lying in a horizontal strip, a contradiction 
with Corollary ( 3.15D (ii) and the local optimality of F. 



Then 9i = 9ii = 0, and, by Remark |4.2| , there exists two points A, B such that the 
segments -S] and {A, Ar} are horizontal, with length greater than or equal to 1, and 

they are contained in F. In order to complete the proof we have to show that |-B,^] is 
a horizontal segment. Assume by contradiction that this is not the case, that is B ^ A', 
where A' is the intersection of the line containing fA, ArJ with Dr-i- In this case the length 
of the polygonal line |Pr-i7^'j^rl is less than the length of any curve joining Br-i with 
Ar and containing [P^-ii^l and ^,.], in contradiction with the local minimality of 
F. □ 

Definition 5.2. An S's-path is a Snell path joining the point (2m + k, k) to the point 
(2m + A; + 3, A; + 1) for some m, /c G N. We shall denote its normalized length by 

As = /(2, /3) = , ^ + , -2-^2, 
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where cjs = a (2, (3) is implicitly defined by 

2(73 , 0-3 



1 = 0. 



^3 ^P' - 4 

The optical length of an S'3-path will be denoted by A3 = A3 + 2 + ^/2. 

Remark 5.3. It is straightforward to check that -v/lO < A3 < |3^/\Q, and, by the very 
definition of A3 = 2 + when [3 = Moreover, l/^/H) < 0-3 < [3/^/lQ. 

The S3-paths will play a fundamental role in the analysis of the geodesies for ^3/2 < 
P < (3q (see Theorem 5.4 below). Namely for /? in this range the 5'3-paths have the minimal 
normalized length among all the Snell paths starting from a light vertex (2m + k, k) and 
reaching a point on the line y = k + 1 (see Corollary |3.15 ). 

Theorem 5.4. Let T he a geodesic as in (fg^j. If \/3/2 < [5 < then the points A^, 
are light vertices for every k = 0, . . . ,n. Moreover Bk-i is connected to either by an 
horizontal segment or by an S^-path for every k = 1, . . . ,n. 

Proof. By Remark 15 we can assume, without loss of generality, that Aq = Bq = (0,0). 
Let us define 

(63) i = min{r G {1, . . . , n}; Ar is a light vertex} . 

Notice that the index i in (|6^) is well defined, since as a consequence of Proposition |4.3| (i) 
and (ii) at least An is a light vertex. Moreover, if i > 2, then F cuts every light diagonal 
Dfc, k = 1, . . . ,i — l, that is, the points A^ and Bf^ coincide and they are not light vertices. 

Let us denote by T' the portion of F joining Aq = Bq = (0,0) to Ai. We are going to 
prove that 

i = 1, and 

^ ^ F' is either an horizontal segment or an S'3-path. 

Once these properties are proved, then Bi is a light vertex and, repeating the procedure 
n times, we reach the conclusion. 

Let m G N be such that Ai = {2i + m, m). If m = 0, then Bq is connected to Ai by an 
horizontal segment, so that i = 1 and ( |64[ ) holds. If m = 1, then by Corollary 3.15| (ii) and 
the local minimality of F, we conclude that i = 1 and F' is and S'3-path. It remains to 
prove that the case m >2 cannot happen. 

For every k = 0, . . . ,i let r^ = 1, . . . , m be such that ^4^ = {2k + r]k,r]k) G F^^., where 
Fr is the Snell path defined in (|6T|). Let us denote by Ck = {xk,rk — 1) and = (x'^j^^) 
respectively the starting and the ending point of F^^. , let 6^ = x'^, — Xk be the thickness of 
Fj.^, and define 

t^ =2k + rk-l-Xk, t^ = x'i,-2k-rk, hk = Vk - [m\ , 

so that 6k = t^ + t^ + 1 (see Figure p^). 

For A; = 0, we have that tq = 1, Cq = Bq = (0, 0) and Cq = (5o, 1). By construction we 
have (5o > 1. We claim that 2 < 5q < 3. Namely, by Corollary 3.15| (v) and (iii), we have 

(Cf^iTi) = l{5o,P) + V2 + So-l>V2 + So-l = £^([0, (1, 1), {60, 1)1) , if <5o G (1, 2], 

|/:;3(Fi) >/(2,/3)+^/2 + 5o-l>/:/3(m(3,l))Ul(3,l),(5o,l)l) , if5o>3, 

hence, by the local minimality of Fi, we cannot have neither 1 < 60 < 2 nor 60 > 3. For 
k = i the same arguments show that 2 < 6i < 3. Hence, from Proposition 4.6, we have 
that 1 = To < vi < ■ ■ ■ < ri = m. 
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Figure 10. 



Finally, by Lemma |3.17 , 

>[(t^,/3,l), l{t+,(3,l-hk)> 1(4,(3,1), k=l,...,i-l, (i>2) 

where I is the function defined in (|55|). Moreover, being F a geodesic, we have l{t'^,(5, h^) < 
and 1(4, P,l — hk) < 0, so that by Corollary |3.15| (v), we have t^,4 > 1- On the 
other hand, from Corollary 3.15 and Proposition |4.6| , we cannot have > 2 or > 2, 
otherwise F^^. would intersect more than one light diagonal. In conclusion, we have that 
t]:,4 e (1,2), so that <5fc G (3,5), A; = l,...,i-l (i > 2). 
Let us define 

fc=0 

By construction and from the estimates above we have 

A<2i + m, m>i + l, 3i + 1 < A < 5i + 1. 
It is straightforward to show that 



cJr'WjTrA >£^(I(A,i + l),(2i + 

\ k=0 ) 



m,m) 



so that 

JOfsir') > E^M^rJ + v/(2i + m-A)2 + (m-i-l)2. 

k=0 

£^(Fj.j.) can be estimated using the function I. For A; = and k = i we have that 

, , Cl3(Tr,) = l{5o- 1,(3,1) +60-I + V2, 

(65) , 

£^(F,J = /(5i - 1, /3, 1) + - 1 + ^ , 

while, for k = 1, . . . ,i — 1, (i >2) 

(66) Cf3iTr,)=l{q,(3,hk) + li4,P,l-hk)+5k-l + V2. 
From Lemma ^!J(iii) and Lemma |3.17l we have that 

lit, (3, h) > kit, (3, l){h - 1) + ~k{2, /?, l){t - 2) + A3 

for every (t, h) E (1, 2) x [0, 1]. (We recall that A3 = l{2, (3, 1), see Definition |5^) . Moreover, 
from (|56| ) and (|5^) we have that 



lt{2, (3, l) = Jl- - 1, lh{t, (3, 1) = VI - a^{t, (3, 1) + d{t, P,l)-V2. 



where (T3 = d{2,(3, 1), see Definition |5.2| . Since t 1— > d{t,(3, 1) is a decreasing function for 
t > 0, and the map s 1— > Vl - s2 + s-^/2 is increasing for < s < 1/ ^/2, we finally obtain 



(67) [(t,/3,/i) > A3+ Vl-^i-l (t-2)+ Jl-ai + a3-^/2 [h - 1 
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for every {t,h) G (1,2) x [0,1]. From (||), (||) and (|67|) we obtain 



Ci3{Tro) > Xs + [ Jl - ai - 1) {5o - 3) + 6o - 1 + V2, 



(68) CpiVrJ > Xs + { Jl - ai - 1] {6^ - 3) + 6^ - 1 + V2^ 



{k = l, 

and 
(69) 











(\/l- 




that 





^ /:/3(r,J > A^l - al + as + (2V2 + 4 - 6,Jl - - as + 2X3 

k=0 



Cf3(r') > AJl-al + a3+ 2V2 + 4-6A/l-a|-a3 + 2A3 i 



+ J{2i + m - A)2 + (m - z - 1)2 . 



Let us now consider the path T" starting from the origin, obtained by the concatenation 
of i 53-paths and the segment connecting the point (3i, i) to Ai = {2i + m,m). Since this 
segment connects two points on the hght diagonal Di, its length is (m — i)V^, so that 

£^(r") = iAs + (m - i)V2 = i{X3 + 2 + V2) + {m- i)V2. 

We are going to show that Cfs{T") < £^(r'), in contradiction with the local minimality of 
r'. We have that 



£^(r') - £/j(r") > {2 + V2-3^l-al-a3 + X3)i- {V2-^l-al-a3 

(70) 



+ V^i + M2 - wyi -f^i -M2 i^V2- yi -^3 

where 

Hi = 2i + m — A > 0, fi2 = 'm — i — 1>0. 
Since 3i + 1 < A and m > i + 1 we have that 

(71) < ^1 < /i2. 
Moreover 

(72) ^fil + - ^1^1 - al - (V2 - Y^l - ai^ =1^2^ 
where 



ifis) = Vl + s^ - s ^l-al -[V2- ^l-a: 



Since < 0-3 < l/\/2, we have that l/\/2 < y 1 — o"! < 1. It can be easily checked that 
ip'{s) < for every s £ [0, 1], hence 

(73) = ip{l) < ^{s), VsG[0,l). 

From (^, (InD, (|72D and (P) we thus get 



(74) Cpir') - Cp{r") > 2 + 72 - 3J1 - ai - + A3 U - V2 - Jl - ai - ^73 . 
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We claim that 

(75) 2 + \/2-3 



0-3 + A3 > \/2 - J 1 - 0-| - (73 > . 



The second inequaUty in (^) easily follows from the fact that < (T3 < 1/ \f2. Concerning 



the first one, by the very definition of A3, and since h ^ }? j 
function in [1, +cxd), we have that 



2 + ^/2-3a/i 



-2s \ 



2 /32 



1 



> -2J1 



+ 



+ 



= -V2 



fT5 



fjg IS an increasing 



0-3 



3/2 



1 



fT5 



3/2 



a.. 



It can be checked that ip is strictly increasing in (0, 1). Since, by Remark 5^, (T3 > l/ylO, 
we have that 

, , , J I \ VW 3^10 



15 2\/l4 



V2>0. 



and (75) is proved. 



imply that CfsiT') 



□ 



Since i > 1, it is straightforward to check that (74) and (^ 
£/3(r") > 0, in contradiction with the local minimality of V . 

Theorem 5.5. Let T be a geodesic from the origin to a light vertex ^ = {x,y), with 
< y < X. 

(i) Iff3> f3S, then CpiV) = x + ( ^2 - l)y. 

(ii) // ^3/2 <P<PS> tf^en 



^/3(r) 



X + (A3 - 3)y, 
A3-V2 , 3V2 ■ 



■ X + 



A3 



if0<y<x/3, 
y, ifx/3<y<x, 



where A3 = A3(/3) is the length of an S^-path, introduced in Definition 5. "4. 
Proof, (i) For /? > /3q it is a straightforward consequence of Theorem |5.l| . 



(ii) Let us consider the case \/3/2 < (i < (3q. From Theorem 5.4 we know that F is the 
concatenation of S'3-paths and segments joining light vertices, lying on light diagonals or 
on horizontal lines. Hence we have that 

£^(F) = tA3 + r + d^/2, 

where t is the number of S's-paths, r is the number of unit horizontal segments, and d 
is the number of diagonals of light squares. It is clear that the three numbers t,r,d G N 
must satisfy the constraints 

3t + r + d = x, t + d = y, 

so that 

£^(F) = (A3 - 2 - V2)t + x + {V2-l)y. 

Since A3 < 2 + \/2, L is minimized by choosing the largest admissible value of t, which is 
y if y < x/3, and (x — y)/2 if y > x/3. (We remark that {x — y)/2 is an integer number, 
since the point P = {x,y) is a light vertex.) In conclusion, if y < x/3 we choose t = y, 
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d = and r = x — 3y, whereas if y > x/3 we choose t = {x — y) /2, d = (3y — x)/2, r = 0, 
obtaining (ii). 

Finally, if /5 = Pq, it is straightforward to check that the lengths of geodesies can be 
computed indifferently as in (i) or in (ii). We remark that, in this case, these two formulas 
give the same result, since A3 = 2 + ^2. □ 

One may wonder if the previous characterization of the geodesies for \/3/2 < (3 < (3q 
remains valid for (5i<(5< (3^. The following example shows that this is not the case. 



Example 5.6. Let V = [(0, 0), (1, 1)] U S'((l, 1), (4, 2)). By Theorem |5j, if < /? < /Jg, 

then r is a geodesic joining the origin to the point £^ = (4,2). Given t G [0,1], let us 
consider the curve 

r(t) = 5((0, 0), (1 + t, 1)) U S{{1 + t, 1), (4, 2)). 

We have that 

L{t) ■.= Cp{T{t))=l{t,0) + l{2-t,(3). 
Moreover, for t = we have r(0) = T, L(0) = A3 + \/2, and 

L'(0) = hm M^iM = /+(o,/3) - ;-(2,/3) 

(76) * 

r 

2 

We recall that, for a given /? > 1, (T3 is the unique zero in (0, 1) of the function 

, ^ 2a a 



/5^-^-Vi--i- 



Moreover, 5 is a strictly monotone increasing function in (0,1), with (^(0) = —1 and 
5(5) — > +00 as s ^ 1~. For /? < ^3/2, let us compute 



Since 



^3_2/?2 V (4/32 -3)3/2 (2/32 _ 1)3/2^ 

the map is strictly monotone decreasing in [1, ^^3/2], with V'(l) = 1 a-nd ^ ^^^3/2^ = — 1. 

The unique zero of in (1, ^3/2) is /3 ~ 1.17868, and /5 > /3f ~ 1.06413. Hence, for 
1 < /3 < we have that 0-3 < v^3/2 - /?2 and, by (|76|), L'(0) < 0. In conclusion, if 
1 < /3 < /3, for t > small enough we have that L(t) < L{0), and F is not a geodesic. 

6. The homogenized metric 



As a direct consequence of Theorem 5.5, we obtain the complete description of the 
homogenized metric for /3 > a/3/2. In the general case we discuss the regularity of the 
homogenized metric. 

In order to make some usefull reductions, we need two remarks on the distance d|(0, ^) 
defined in (^). 

Remark 6.1. Since |C — ??| < (i|(??, ^) < (3\^ — ri\, it can be easily seen that 
(77) <Pp{0 = ^lim d^(%, e^) , V C G M2, V6 ^ e , Vr?, ^ . 
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Figure 11. The homogenized unit ball for /3 > (left) and \/3/2 < (3 < 
PS (right) 



;§,§). From r 



we have that 

l>2 



Remark 6.2. For every e > let ??e 

Since the map ^ i— > ^+r/e) is symmetric w.r.t. the coordinated axis and the diagonals 

passing through the origin, it is clear that has the same symmetries. 

In what follows, given (x, y) G we denote 

M = max{|j;|, |y|}, m = min{|x|, |y|}. 

Theorem 6.3. For every (x, y) G M? the following hold. 

(i) IfP> then ^^(x, y)=M+{V2- l)m. 

(ii) // <P<Po> *^en 

'm + (A3 -3)m, 
^/3(2;,y) = <( A3 



, 3^- A3 
iVi H m, 



if 3m < M, 
if3m>M, 



2 2 

where A3 = A3(/3) is t/ie length of an S^-path, introduced in Definition 5.L 

Proof. By Remark |6.2| it is enough to consider the case < y < x, so that M = x and 
m = y. Let = (x, y) and, for e > 0, let us define 



3 



n 



x-je 
2e 



= {xs,ye) = {(2n + j)e,j£). 



Then 

\y - 2/el < \X- Xe\ < 2£, 

SO that 1^ — < £^/5. Moreover (i|(0, ^e) is explicitly computed in Theorem 5.5. Since 
^, by (|77|), the conclusion follows. □ 



In the general case we have the following result. 

Theorem 6.4. Let P > 1 be given, and let kc{P) be the number defined in Definition^ 
Then <P(}{x,y) = M + {l{2kc{P), p) + V2 - l)m for every (x, y) G belonging to one of 
the cones {(2/cc(/3) + l)|y| < |x|} or {{2kciP) + < |y|}- 



Proof. From Remark 3^ it is enough to consider the case < {2kc{P) + l)y < x. Moreover, 
we can assume that y > 0, the case y = being trivial. Since the homogenized metric 
depends continuously on /?, it is not restrictive to assume /3 7^ in such a way that 

K2kciP),P) < l{t, P), V t > 0, t / 2k,{p) 
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(see Corollary 3.15| (iii)). For every < e < y, let = {xe,ye) be the nearest light vertex 



to ^ = (x,y) below the line x = {2kc{P) + l)y. Then = ((2n + for some j > 1 

and n > kc{P)j, and — > ^ as e tends to 0. 
We claim that 

(78) 4(0,e.) = x, + (/(2/c,(/3),/3) + \/2-l)y,, y < e < y 



so that the result will follows from (|77|). In order to prove the claim, after a scaling, we 
have to depict a geodesic F joining the origin to the point (2n + in the standard 
chessboard structure. Let us define the class S of all Snell paths joining the light vertices 
(2m + r — 1, r — 1) and (2m + 2kc{(3) + r,r), m, r G Z. We are going to show that F has to 
be the concatenation of j Snell paths in S and of horizontal segments. As a consequence, 
since the length of any path in S equals to l{2kc{P), P) + 2kc{f3) + \/2, whereas the total 
length of the horizontal segments is 2(n — kc{P)j), we obtain that ( [78|) holds. 
Let us define the paths Fi, . . . , Fj as in 

(13) • For every r = 1, . . . , j let us denote by 
{xr,r — 1) and (x^,r) the endpoints of F^, and let Thanks to Lemma |3.16| , 

we have that 

£^{Tr) > l{2k,{(]),P) +Tr-1 + V2, 

for every r = 1, . . . ,j. Then we get 

j j 
Cf3iT) = J2 ^fsi^r) + 2n + j-J2^r>2n + j + (/(2fc,(/?), /?) + V2 - l)j 

r=l r=l 



and again by Lemma 3.16 the equality holds if and only if = 2kc{P) and F,. G S for 



every r = 1, . . . ,j. □ 

Corollary 6.5. For every (3 > 1 the unit hall of the homogenized metric is not strictly 
convex, and its boundary is not differentiahle. 

Remark 6.6. The presence of faces in the optical ball corresponds to nonuniqueness of the 
geodesies. More precisely, if F is a face of positive length, and C = {Ary; r] G F, A > 0} is 
the corresponding cone, then for every ^ G C, a function u £ AC{[0, with u{0) = 0, 
u(l) = ^, parameterizes a geodesic if and only if u'{t) £ C for a.e. t G [0, 1]. 

Namely, there exists p G such that ^^{rj) = {p, rf) for every rj G C, and $/3(r/) > {p, rj) 
for every r/ G \ C. Hence, if u'{t) G C for a.e. t G [0, 1], then we get 

Cf^{u)= j\^{u'{t))dt= f\p, u'it))dt = <^^iO, 
Jo Jo 

whereas C^°^{u) > ^/siS,) whenever {t G [0, 1]; u'{t) C} has positive measure. 

As a final remark, let us consider the chessboard structure corresponding to the upper 
semicontinuous function 



(79) dpiO 



p ifee ([0,1] X [1,2]) U ([1,2] X [0,1]) 
1 otherwise 



which differs from the standard chessboard structure defined in ([T|) by the fact that d = {3 
instead of 1 on the sides of the squares. 

In this way we obtain a new family of length functionals In this case the existence 
of a geodesic joining the origin with a point ^ G is not guaranteed. For example, if 
i = (e, 0), we have Z%{u) > e for every u G AC{[0, 1], M^) such that u{0) = and u(l) = 
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On the other hand, we can construct a minimizing sequence («n)n such that Cp{un) s 
for n — +00, defining 




{st,2t/n), if [0,1/2], 

{et,2{l-t)/n), if [1/2,1]. 



Nevertheless, the F-Iimit with respect to the L^-topoIogy of the functionals {Cp) coincides 
with the F-Hmit of the functionals Namely, the hminf inequality is certainly 

satisfied since £| > On the other hand, given u G ^C([0, 1],M^) and a recovering 
sequence (ug) for we can construct a recovering sequence (tie) for (>C^) in the following 
way: for a given e, we obtain modifying Us in the region where Uelt) belongs to the 
set S of the sides of squares, in such a way that J0.p{ue) < jC-p{ue) + £, and the set 
{t G [0, 1]; Ue G S} has vanishing Lebesgue measure. 
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